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The effect of geometrical frustration in a two-dimensional 1 /4-filled strongly correlated electron 
system is studied theoretically, motivated by layered organic molecular crystals. An extended Hub- 
bard model on the square lattice is considered, with competing nearest neighbor Coulomb inter- 
action, V , and that of next-nearest neighbor along one of the diagonals, V' , which favor different 
charge ordered states. Based on exact diagonalization calculations, we find a metallic phase sta- 
bilized over a broad window at V' ~ V even for large Coulomb repulsion strengths as a result of 
frustrating the charge ordered states. Slightly modifying the lattice geometry relevant to the actual 
organic compounds does not alter the results, suggesting that this 'quantum melting' of charge order 
is a robust feature of frustrated strongly correlated 1/4- filled systems. 
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Effects of geometrical frustration in strongly correlated 
systems have been studied actively, where novel exotic 
states may arise as a result of competition between differ- 
ent 'conventional' ordered phases. Such possibilities have 
been explored extensively in spin systems, for example, in 
spin- 1/2 models with frustrated antiferromagnetic (AF) 
exchange interactions of nearest neighbor (NN) and of 
next-nearest neighbor (NNN)A There, conventional Neel 
states can be destabilized due to the competition between 
different AF patterns and quantum fluctuations lead to 
exotic states such as spin liquids. In contrast, works de- 
voted to the charge degrees of freedom are limited and 
little is known on such situations. One can consider an 
analogy between the classical spin systems, i.e., Ising sys- 
tems, and the charge ordering (CO) systems: 2 " However, 
when CO melts due to quantum fluctuation, nature of 
this transition will be different from spin systems because 
of Fermi statistics. Furthermore, quantum melting of CO 
is of particular interest since CO states have now been 
found in various materials such as transition metal ox- 
ides^ and organic conductors^ including systems suscep- 
tible to frustration such as in a triangular lattice system 
Naa;Co02 recently attracting interest £ 

Our motivation here is CO systems having frustrated 
lattice structures, namely, quasi two-dimensional (2D) 
organic conductors 0- and a-ET 2 X (ET = BEDT-TTF, 
X: monovalent anion). These systems have recently 
been studied extensively, displaying a subtle competition 
between CO, metallic and superconducting phases 
Their electronic states have been studied theoretically 
by using 2D extended Hubbard models which include 
Coulomb repulsions of not only on-site, U, but also in- 
tersite, Vij. The latter plays a crucial role for CO at 



3/4-nlling (or equivalently 1/4- filling of holes) AiS So far, 
this kind of models have been studied for simple lat- 
tice structures such as one-dimensional (ID) chain^ lad- 
der^ 2D square lattice ; 14 ! 15 and so on. In all of these 
works, Wigner crystal-like CO states become stabilized, 
in general, due to the effect of Vij between NN sites when 
U is large. 

Recently, however, there appear experiments in these 
organic materials which are difficult to interpret within 
naive understandings of such CO phase transitions. In 6- 
ET2RbZn(SCN)4, when CO was quenched by rapid cool- 
ing, a glassy CO state is foundAi& Even above the tran- 
sition temperature where the system shows metallic re- 
sistivity^ the dielectric constant is unexpectedly found 
to have a frequency-dependence characteristic of insula- 
tors in members of 6-{FIT)2X ^ while an extremely slow 
dynamics of CO is found in 13 C-NMR in the above 0- 
type compound as well as in cv-ET 2 l3ii& These suggest 
effects of geometrical frustration in the charge degree of 
freedom, which we investigate in this paper. 

Fortunately, the modelling of organic systems is very 
simple and reliable. The transfer integrals between 
the molecules, ty, can be obtained from the extended 
Hiickel method^i In the 9- and cv-type structures, the 
ET molecules are arranged in 2D planes which are close 
to triangular lattices. However, the network of t%j can 
be approximately described by a square latticed On the 
other hand, the values of Vij are comparable on every 
bond of this triangular lattice^ which introduce frustra- 
tion in this system. In the previous calculations on 2D 
extended Hubbard models for these materials, 10 U and 
Vij were treated at the mean-field level and only insulat- 
ing phases were found for large interactions. However, 
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we expect that a 'quantum melting' of CO takes place 
due to frustration, and that a metallic region is stabi- 
lized even in the strongly correlated regime. Actually, in 
a frustrated ID model, destabilization of CO state due to 
quantum fluctuation was shown numerically. 19 Therefore 
it is of fundamental interest to proceed beyond mean- field 
in those 2D models. 

To study this problem, we first consider a simplified 
model for the 2D organic systems based on the consider- 
ations above, i.e., the 1/4-filled extended Hubbard model 
on a square lattice with NN and NNN Coulomb repul- 
sions: 

ff = E (4v C fc + c]a C i°) + U ^T^l 

(ij) (ay 

Here, c ia and (cj CT ) creates (eliminates) an electron of 
spin a at site i, n; L(J — c? ia Ci a and rii = n,-f + n,-^. t and V 
are the transfer integrals and the NN Coulomb repulsion 
along the bonds of the square lattice (ij), respectively, 
and V is the NNN Coulomb repulsion along one of the 
diagonals (ij)'- 

Some insight of the CO states in the above model can 
be attained by considering the following two limits for 
U/t 3> 1, where doubly-occupied sites are suppressed: 

Case A. V/V 3> 1: An insulating state due to a 
checkerboard pattern of CO with alternating charge in 
every other site (wave vector q = (n, 7r)) is realized. The 
effective model for the spin degrees of freedom on the sites 
occupied by an electron is a 2D spin-1/2 AF Heisenberg 
model with NN J and NNN J' along both the diagonals. 
The estimated exchange couplings from fourth order per- 
turbation in t show J> J' and then q = (ir/2, ir/2) spin 
ordering is expectedii^ 

Case B. V /V 3> 1: The system is insulating but 
with q = (0, 7r) CO forming rows or columns of occu- 
pied sites alternating with empty ones in a stripe-type 
manner. The spin degrees of freedom is described by a 
quasi- ID AF Heisenberg model with J along the chains 
coupled through the interchain interaction J'. Therefore 
the AF spin ordering of q = (it, tt/2) is expected. 

When V'/V ~ 1, the frustration is expected to play a 
major role. Actually, if t — and U is large enough, the 
model is equivalent to the 2D Ising model with couplings 
V/2 and V'/2 as can be easily seen from Eq. (QJ. In this 
case the ground state for exactly V — V is always gapless 
which corresponds to the metallic state in the present 
model. Below we will analyze Eq. Q by Lanczos exact 
diagonalization on an L — 16 cluster to take account of 
quantum fluctuation. 

In order to determine whether the system is insulating 
or not we have calculated the Drude weight^: 
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FIG. 1: Ground state phase diagram of the 1/4-filling ex- 
tended Hubbard model Eq. Q for U = Wt. The boundaries 
are extracted from exact diagionalization calculations of the 
Drude weight on L — 16 site cluster (see Fig 01 - The label M 
stands for metallic, and COI for a charge ordered insulator. 



where Eq and E n are the ground state and excited 
state energies of the system, respectively. T is the ki- 
netic energy operator (first term in Eq. l(Tjl). j x the 
current operator in the x-direction and e the electron 
charge. To explore the electronic states we have also cal- 
culated the charge and spin structure factors, C(q) = 
i E . . e iq.R y {n . nj) and s (q ) = i £. . e i*y ( S ?S?), re- 

spectively, where is the vector connecting two differ- 
ent sites and Sf is the z-component of the spin operator. 

The ground state phase diagram on the V — V plane 
for a fixed U — lOt is shown in Fig. 2] which includes CO 
insulating phases and a metallic phase. At weak coupling 
with small (V/t, V'/t) a metallic phase is stabilized as ex- 
pected from the large fraction of charge carriers available 
in the system. Fixing V' = and increasing V induces 
a metal-insulator transition at V ~ 2t to the checker- 
board CO phase (Case A discussed above)^ On the 
other hand, if we fix V = and vary V' , a metal-insulator 
transition to the stripe CO phase occurs at V' ~ At (Case 
B). Interestingly, the metallic phase appearing at weak 
coupling extends out along the V = V line up to strong 
coupling in a robust manner. 

In Fig.|2the Drude weight for different values of V and 
in Fig. [3] the structure factors for V = 3t as a function 
of V', both for fixed U = lOt, are plotted. Fixing V to a 
large value (V > 3t), we find that the checkerboard CO 
state with AF ordering is stabilized for small values of V, 
as observed from a negligibly small Drude weight and the 
large values of C(ir,ir) and S(ir /2,ir/2). This insulating 
phase is destabilized with increasing V and a transition 
from the CO insulating phase to a metallic phase oc- 
curs as can be seen from the rapid increase of the Drude 
weight with V' . In this region, C(tt,tt) and S(ir/2,n/2) 
are suppressed with V' and a crossover to the stripe CO 
state is seen, characterized by the large value of C(0,7r) 
and S(n, tt/2) at large V. When V - V, C(q) becomes 
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FIG. 2: The Drude weight calculated by exact diagonalization 
on L — 16 for the 1/4- filled extended Hubbard model Eq. Q 
for U = Wt and different values of V. It is rapidly enhanced 
as V' is increased with V fixed signalling the presence of a 
metallic phase. Inset shows for V = 3t the variation under 
including a diagonal hopping, t' = O.lt, to the model, which 
is small. 



featureless indicating that neither the checkerboard nor 
the stripe CO state is preferred, as expected. 

The existence of the intermediate metallic phase found 
here should be robust enough in the thermodynamic 
limit, although the cluster size investigated is limited to 
L = 16. 21 Whether the CO quantum phase transition is 
always accompanied with a metal-insulator transition or 
not is not clear from our results. As a matter of fact, 
sudden changes of S(tv/2, 7t/2) and the Drude weight at 
V'/t = 1 indicate a level crossing to a metallic state. 
However, C(tt,tt) still dominates over the other q val- 
ues suggesting the presence of CO metallic states near 
the boundary. Calculations on larger systems and finite 
size scaling are neccesary to determine the CO transition 
lines more accurately. 

This frustration induced metallic phase is found to be 
robust against introducing a small hopping amplitude t' 
along the diagonal direction (ij)', relevant to the actual 
0- and a-(ET) 2 X materials^ In the insets of Figs. H 
and|3 comparison between the square lattice case Eq. (JTJ 
and the results for the model where t' — O.lt is added are 
shown. The metallic phase is stabilized in a similar V 
range, and C(7r,7r) (C(0, 7r)) diplays similar suppression 
(enhancement) as V' is increased. 

Finally, we discuss the evolution of the optical spec- 
tra^ upon these crossovers, as shown in Fig. 0] for fixed 
U = lOt and V = 3t for different values of V. In the case 
of V' — (top frame), the system is in the checkerboard 
CO insulating phase, and a gap in the spectra with a peak 
position centered around 8t is observed. This is naturally 
expected since moving a charge from the checkerboard 
CO state costs energy of 3V (= 9t in this case) in the 
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FIG. 3: The charge and spin structure factors, C(q) and S(q), 
for U = lOt and V = 3t as a function of V' for wave vec- 
tors relevant for the two CO states discussed in the text. A 
crossover corresponding to the Drude weight is seen. Inset 
shows the variation of C(q) when small diagonal t' = O.lt is 
included in the model, where the change is small. 



atomic limit On the other hand, in the stripe CO in- 
sulating phase, there are two possible ways of moving a 
charge, which are perpendicular and along the stripes, 
with the estimated energy costs of 2V' — V and U — V, 
respectively. These lead to a rather broad incoherent 
band around 7t in the bottom frame of Fig. 0] in the 
case of V = 5t, as expected. Apart from this, sharp 
peaks at about At appear, which can be attributed to 
the ID nature of the stripes. Indeed, excitonic features 
appear inside the Mott gap in the optical spectra of a 
1 /2-fflled ID extended Hubbard chain for V < U/2^ 
which should be an effective model for each stripe. 

In the metallic phase at intermediate V'/V (see middle 
frame of Fig. |3J, a redistribution of the spectral weight 
occurs. Spectral weight associated with incoherent tran- 
sitions due to U and Vij is transfered to lower energies 
giving rise to the Drude peak and a mid-infrared band 
together with a low energy feature appearing at its lower 
edge. The overall appearance of the optical conductiv- 
ity is similar to the one observed for a 1/4-filled metal 
close to the checkerboard CO transition^ induced by V 
(V' = 0). We stress again that this is nontrivial since the 
interaction strengh is much larger here in the frustration 
induced case. However, interestingly, the mid-infrared 
band reaches higher energies making the spectrum some- 
what broader while the feature at its lower edge appears 
at larger energies. All these 'exotic' features appearing 
in the frustrated metallic phase can be attributed to the 
dynamics of quasiparticles inmersed in a charge fluctu- 
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FIG. 4: Evolution of optical conductivity ct{u>) of the model 
Eq. Q as the system is driven across the frustration induced 
metallic phase from the checkerboard to the stripe CO phases. 
Spectral weight is transfered from high to low energies as the 
system is driven through the metallic phase forming a mid- 
infrared band with a feature appearing at its lower band edge 
and a strong Drude component. A Lorentzian broadening of 
t; = O.li is used to smoothen out the delta peaks. 

ating background and would be missed in a mean-field 
analysis. 

In the actual organic compounds, 9- and a-(ET)2-X", 
U ~ 1 eV and t ~ 0.1 eV which correspond to U ~ lOt, 
and V' ~ V with the value of ~ U/3 is realized^ there- 
fore the system is inherently located in or near to the 
frustration induced metallic phase. In reality, however, 
the more lower symmetry for ty and/or the coupling to 
the lattice degree of freedom, which we do not include in 
our model, would push the system toward CO states. As 
a matter of fact, large displacements of the molecules are 
observed in ev-(ET)2l3 at the CO transition temperature, 
and in 0-(ET)2Y a structural change lowering the crys- 
tal symmetry takes place concommitantly with the CO 
phase transition^ One may consider these experimen- 
tal facts as suggesting that the CO phase is realized by 
relaxing the frustration effect seen in our calculations, 
since the interaction strength is already strong enough 



to drive the CO state. Another aspect worth considering 
when comparing our results to experimental data is that 
the CO states actually observed in ET compounds have a 
different charge pattern as the ones discussed here^ cor- 
responding to a 'zigzag'-type pattern in our model. We 
mention that this phase is also found in our exact diag- 
onalization study, e.g., for appreciably large V', larger 
than 6t for U = lOt and V = 3t. We also note that the 
CO pattern corresponding to the checherboard pattern 
has been observed in a similar material with the 0-type 
structure, i.e., (9-(BDT-TTP) 2 Cu(NCS) 2 . 23 The precise 
determination of the phase diagram including such dif- 
ferent CO states in comparison with each material, as 
well as more detailed investigation of the frustration in- 
duced metallic phase are left for future studies. 

In conclusion, we have found a metallic phase induced 
by frustration in 2D 1/4-filled systems which is robust 
against changing the lattice geometry from a square to 
an anisotropic triangular lattice relevant to layered or- 
ganic conductors. Our phase diagram is similar to that 
of the extended Hubbard model on the ID zigzag chain 
structure^ as well as to that of a frustrated 2D Bose- 
Hubbard Hamiltonian where a superfluid phase is stabi- 
lized between two charge modulated phases^ From our 
findings, it is natural to conjecture that 'quantum melt- 
ing' of CO states due to geometrical frustration is rather 
general and the metallic phase near the CO phase in the 
actual compounds with geometrically frustrated struc- 
ture is relevant to such situation. An interesting issue to 
address in future is whether superconductivity can exist 
between the frustration induced metallic phase and the 
CO phase mediated by charge fluctuations appearing at 
the phase boundary. This kind of possible superconduc- 
tivity has been studied in connection to organic conduc- 
tors^ and to N& x Co0 2 & 
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